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Abstract 

We study the phase structure of a phantom tethered surface model shedding hght on the internal 
degrees of freedom (IDOF), which correspond to the three-dimensional rod like structure of the lipid 
molecules. The so-called tilt order is assumed as IDOF on the surface model. The model is defined 
by combining the conventional spherical surface model and the XY model, which describes not 
only the interaction between lipids but also the interaction between the lipids and the surface. The 
interaction strength between IDOF and the surface varies depending on the interaction strength 
between the variables of IDOF. We know that the model without IDOF undergoes a first-order 
transition of surface fluctuations and a first-order collapsing transition. We observe in this paper 
that the order of the surface fluctuation transition changes from first-order to second-order and to 
higher-order with increasing strength of the interaction between IDOF variables. On the contrary, 
the order of collapsing transition remains first-order and is not influenced by the presence of IDOF. 
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I. INTRODUCTION 



The crumpling transition has long been interested in membrane physics and in bio 



physics 
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The curvature model of Helfrich, Polyakov and Kleinert 
membranes was found to undergo first-order transitions on spherical and fixed connectivity 
surfaces by Monte Carlo (MC) simulations 

Internal degree of freedom (IDOF) corresponding to the three-dimensional rod like struc- 
ture and electrostatic structures such as a dipole moment of molecules are out of considera- 
tion in the curvature models. The surface models are those defined only by two-dimensional 



differential geometric notions 111]. For this reason, the thermodynamic properties of the 
models can easily be accessed so far in theoretical/numerical studies 
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However, three-dimensional structure of molecules is considered to play important roles 
in specific phenomena in membranes. On the Langmuir monolayer, photoinduced traveling 



waves were observed experimentally 



17l |. The traveling wave is carried by the rotation of 



the molecular azimuth, where the tilt angle is kept constant. This clearly indicates that the 
molecules tend to align to each other. The chirality of membranes is also considered to be 
connected to the tilt of lipids. The directional order- disorder transition corresponding to 
such three-dimensional structure of molecules is the so-called gel-liquid crystal transition, 
which can be observed in bilayers including biological membranes. Lippling transition |l^ 
is also considered to be connected to the IDOF such as directional order-disorder of the 
lipid molecules. Moreover, it is also quite well known that a variety of shapes and topology 
of membranes are both closely related to internal molecular structures. In facts, lamellar, 
hexagonal, and vesicles are understood to be originated from the difference in the shape of 
lipids. 

The tilted molecules and its relation to the shape and the chirality of membranes have 



long been studied 



19|, 



21 



22|. An interaction between tilt order and surface was also 



studied in a model of membrane 
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2J]. An interaction between the shape of membranes 



and the tilt order was studied with the renormalization group strategy. It was reported how 
thermal fluctuations of membranes associate with the strength of the interaction. However, 
a relation between the crumpling transition and the tilt order is still remained to be studied. 

Therefore, it is interesting to see how the crumpling transition depends on such IDOF. 
In order to see influences of the IDOF on the transition, we assume a three-dimensional 
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FIG. 1: (a) Three-dimensional rod like structure of molecules, which is called a director, (b) the 
projected component py of a director p, and (c) the unit vectors m on the triangles. The vector 
m is defined by m = pjj/|pjj |. 

director pj on the triangle i, which is an element constructing a surface. Although IDOF of 
lipids cannot always simply be expressed by the tilt, we simply consider the lipid as a three 
dimensional vector. The directors are drawn schematically as three-dimensional vectors in 
Fig. [T](a). The Hamiltonian corresponding to such IDOF is described by a conventional 
local spin-interaction between unit vectors rrij, which are defined by using the projected 
component of Pi\\ parallel to the triangle i such that mj = pj||/|pj|| | as shown in Fig. [T](b). 
Figure [D^c) shows the vectors rrij on triangles, which are elements of the surface. The 
interaction is identical to the one of XY model if the surface is flat, however, it becomes a 
three-dimensional one on curved surfaces in the sense that the normal perpendicular to the 
unit circle 5*^ (the phase space) varies in R^. 



Our model is similar to the model in 



23l . |24| |. and therefore the interaction between the 



surface and the tilt order is taken into account. It is very interesting that the IDOF of 
molecules interacts with the external degrees of freedom that are the shape of surfaces. In 
many statistical systems the external degrees of freedom and the IDOF are treated inde- 
pendently. We know that spin models such as XY model and Potts model on fluctuating 
surfaces were extensively studied {25!. I26, 2?! . 28 . 29|. however, their IDOF interact only with 



the intrinsic geometry and hence seems to be independent of the shape of surfaces. 

We expect that the model in this paper reveals a non-trivial influence of the tilt or- 
der on the crumpling transition. If it were not for the projection of Pi on the triangle 
i, the variable pj has no connection to the extrinsic geometry of surfaces. However, the 
projected variable rrij obviously interacts with the surface, therefore the shape of surface 
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is influenced by correlations between irij. Conversely, we can also expect that the surface 
fluctuation nontrivially influences the interaction between mj. The so-called KT transition 
of XY model on periodic flat surfaces can be changed into some other transitions or can 

disappear by three-dimensional effects caused by the fluctuation of surfaces or by effects of 
the surface topology. However, in this paper we concentrate on how the tilt order influences 
the crumpling transition in a broad range of interaction strength between the variables mj. 

II. MODEL 

By dividing every edge of the icosahedron into L pieces of uniform length, we have a 
triangulated surface of size = 10L^-|-2 (= the total number of vertices). The starting 
conflgurations are thus characterized by = 12 and Nq = N — 12, where A^^ is the total 
number of vertices with the co-ordination number q. 

The surface model is defined by the partition function 

„/ N 

Z^Y. \{dXieM-S{X,T,xn)], (1) 

m i=l 

S{X, T, m) = 5*1 4- 65*2 + a^a, 

where h is the bending rigidity, a is the coefficient of the XY model, and /' denotes that 
the center of the surface is fixed. S{X, T , m) denotes that the Hamiltonian S depends on 
the position variables X of the vertices, the triangulation T, which is fixed, and the variable 
m. Ylim denotes the summation of the IDOF corresponding to the Hamiltonian ^'3 for the 
XY model, which is defined by 

53 = 5](l-m,-m,), (2) 

where '^3 o^®^ bonds (ij), which are edges of the triangles i and j. The 

vector nij is defined on the triangle i and a three-dimensional unit vector parallel to the 
triangle. As mentioned in the introduction, irij corresponds to the projected component pj y 
of the director p^; mi oc pj y and has values in the unit circle on the plane parallel to 
the triangle i. The director pj tilts a constant angle from the normal of the surface and 
rotates around the normal vector and interacts with the nearest neighbors. This induces the 
interaction described by ^'3. Note also that depends on the curvature of the surface. 
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We emphasize that the model of Eq.([2]) is not identical to the naive XY model, whose 
variables have values in the unit circle 5*^ in R^. Nevertheless, we call the model defined 
by 5*2 in Eq.([2]) as the XY model as mentioned in the Introduction. The reason why we 
call the model as the XY model is because the interaction is almost two-dimensional on a 
smooth two-dimensional sphere, which is locally flat. 

The Gaussian bond potential Si, the bending energy term 5*2 are deflned by 

S^ = J2iX^-X,)^ 52 = 5^(l-n,-n,), (3) 

where ^(j^) in Si is the sum over bond (ij) connecting the vertices i and j, and 

5*2 is also the sum over bond (ij), which is the common edge of the triangles i and j. The 

symbol rij in Eq. denotes a unit normal vector of the triangle i. 

III. MONTE CARLO TECHNIQUE 

The vertices X are shifted so that X' = X+6X, where 6X is randomly chosen in a small 
sphere. The new position X' is accepted with the probability Min[l, exp(— AS")], where 
AS = S'(new) — ^(old). The value of Hamiltonian S^ also changes due to the shift of Xj, 
since the normals of the triangles touching the vertex i change. Then, the new vector is 
obtained by flrstly projecting the old nij to the new triangle and secondly normalizing the 
projected vector to the unit length. The radius of the small sphere for 6X is chosen so that 
the rate of acceptance for X is about 50%. We introduce the lower bound 1 x 10^® for the 
area of triangles; however, no triangle appears whose area is less than 1 x 10~®. Therefore, 
we can say that no lower bound is imposed on the area of triangles. No lower bound is 
also imposed on the bond length. We call a sequential updates of X as one Monte Carlo 
sweep (MCS). 

The vector irij has its value in a unit circle, whose normal is parallel to rij a normal vector 
of the triangle i. The new vector m- is randomly chosen in the circle; is independent of 
the old rrij. As a consequence, we have about 60% ~ 70% acceptance rate for the random 
shift of rrij. The variable rrij can be updated even when X is updated, because the updates 
of X change the normal vectors of triangles and hence rrij. 

Convergence speed of MC simulations for the variables nij is very fast compared to that 
for X, because the phase space of rrij is compact (a circle ^^(c R^) ) whereas that of X 
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is noncompact (R^). Therefore, the update of nij is performed at every 10^ MCS in this 
paper. We consider that nij should be updated after the surface shape is deformed to some 
extent. 

We use surfaces of size A^ = 2562, A^ = 4812, A^ = 8442 and A^ = 14442. The thermahzation 
MCS is sufficiently large: 5 x 10^ ~ 4 x 10^, which depends on the size of surfaces and on 
the values of the parameters b and a. 

IV. RESULTS 

A. XY model on rigid spheres 

In this subsection we present the results of MC simulations for the XY model on rigid 
spheres, in order to see the dependence of the correlation of spin variable nij on the coupling 
constant a. The variable X of the surface is frozen in the MC simulations on the rigid sphere. 

The XY model is defined by the partition function 



on the rigid sphere. The radius of the sphere can be chosen arbitrarily, because the partition 
function Z in Eq.(jlj) depends only on the size N and is independent of the radius of the 
sphere. 

The XY model undergoes the KT transition on a fiat lattice. However, we immediately 
understand that the long range directional order disappears on a sphere, although the di- 
rectional order remains locally on the surface. As a consequence, the KT transition of the 
model is expected to be softened on the rigid sphere. 

In order to understand the configuration at a oo, which corresponds to the zero 
temperature because of the unit of a, we show in Figs. [2](a) and Et^b) two configurations 
typical to such sufficiently large a on the surface of size A^ = 812. The configuration in Fig. 
[2][a) obtained at a = 50 has a pair of vortices on the sphere. Flows of m, denoted by cones, 
emerge from (go into) one vortex and go into (emerge from) the other vortex which is in 
the opposite side of the sphere. Two singular points do not disappear even at sufficiently 
large a because of the topological reason, in contrast to the case of periodic planar lattices. 
The configuration in Fig. EJ^b) obtained at a = 200 corresponds to the low temperature 




(4) 



m 
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FIG. 2: (Online color) (a) The KT type configuration at a = 50, and (b) the low temperature 
type configuration at a = 200, where small cones denote the vectors m. A pair of singular points 
appears on both of the spheres; one of them can be seen in the figures, and the other is in the 
opposite side of the spheres. The surface size is A^ = 812, and the magnetization M/Nt of Eq.(l6]) 
is M/iVT=i 0.57 in (a) and M/iVT~ 0.77 in (b). 

configuration. Two configurations are easily obtained by increasing a from small a sucli as 
a = 5 step by step in tlie simulations. The configuration of Fig. [21(b) is almost stable, while 
that of Fig. [2]^a) is unstable and changes to/from that of Fig. [2]^b). 

We consider that the singular points of m on the sphere have no influence on the KT 
transition. In fact, the low temperature configuration has no vortex as we see in Fig. [2](b), 
and XY model on the sphere is reported to have KT transition although the variables m 
have values in the circle S^{c. R^) in contrast to the model in this paper 25||. We should 
also comment on the problem of frustration, which in general appears to influence the phase 
structure of models defined on triangular lattices. However, the model in this paper is 
defined on the lattice that allows no frustrated configuration. In fact, the variables m are 
defined on the plaquettes (= the faces of triangles), and therefore the interaction between 
m forms the hexagonal (or pentagonal) lattices, which are the so-called the dual lattice. 
Therefore, we expect that the results are not influenced by the frustration. 

Figure [3](a) shows S^/Nb against a over the range of couplings spanning the critical 
region of the KT type transition. The symbols Nsi^'iN — Q) and A^'r(=2A^ — 4) denote the 
total number of bonds and the total number of triangles, respectively. Figure [3](b) is the 
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FIG. 3: (a) The internal energy S^/Nb vs. a, (b) the specific heat Cs-j, vs. a, and the magnetization 
M/Njn vs. a. The model is the pure XY model defined by Eq.([3]) on rigid spheres. Cs-j has the 
peak at a ~ 1.45. Nb and Nt denote the total number of bonds and the total number of triangles, 
respectively. 

specific heat C53 for S3, which is defined by 

Cs-. = ^{{S3-{S3yf). (5) 

The size of surfaces is assumed as = 812 ~ = 4842 in the simulations on the rigid 
spheres. 

We find in Fig. [3](b) that C53 has a peak at a ~ 1.45, where the peak value remains almost 
constant as increases. This implies that the model undergoes a higher order transition or 
the KT transition just like the XY model on the fiat regular lattice. The point to note is 
that lUi variables become relatively ordered at a > 1.45 and relatively disordered at a < 1.45 
on the rigid spheres. 

The magnetization M is defined by 

M=||M||, M = J^mi, (6) 

i 

where is the sum over the triangle i. Figure [3]J^c) shows M/Nt versus a. 

As we have seen in the low temperature configuration in Figj2]^b), M is expected to be 
large at sufficiently large a even on the rigid sphere. For this reason, we expect that M 
reflects ordering/disordering of m although m is defined parallel to the spherical surface. 
On fluctuating surfaces including those in the collapsed phase, M is expected to have a role 
of order parameter of the transition. 
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The variance xm oi M can also be defined by xm = ( {M- {M)Y) /Nt- The KT transition 
in the XY model on the fiat lattice is the one that is characterized by the divergence of xm 
at the transition point. It is interesting to see whether the KT transition persists in the 
model of this paper just as in the XY model in 25|]. However, we do not go into this point 
further as mentioned in the last of the introduction. 



B. Collapsing transition 

In the following we concentrate on the influence of IDOF on the phase structure of the 
surface model. In order to do this, the coupling constant a is fixed to a = 2, a = 1, and 
= 0.5. Since the transition point is a ~ 1.45 on rigid spheres, we expect that the value of 
a = 2 is sufficiently large for the m variables to tend to align themselves even on fluctuating 
spheres. On the contrary, the values a = l and a = 0.5 are considered to be sufficiently small 
to disorder the m variables. 

Before analyzing the collapsing transition, we show some of the quantities such as the 
magnetization M/Nt and the internal energy S^/Nb to get information on the configurations 
corresponding to the m variables. Figures |l](a) and|l](b) show S^/Nb and M/Nj- versus the 
bending rigidity b obtained under the conditions a = 2, a = 1, and a = 0.5. The data 
are obtained on fluctuating spheres. The range of b in each a is the region of collapsing 
transition point be, where the surface collapses at b<bc and swells at b>bc. 

The value S^/Nb — 0.37 at a = 2 in FigJl](a) indicates that the interaction strength a = 2 
corresponds to the ordered phase of the XY model on rigid sphere, which was shown in 
Figsl3](a) and[3]^b). We see also in FiglD^a) that the vectors rrij are almost decorrelated to 
each other at a = 1 and a = 0.5. We get the same information on the configurations of 
rrij from M/Nt in Figj4](c) by comparing the values of M/Nt to the corresponding ones in 
FigHc). 

We find also from FigsJD^a) and|l](b) that both M/Nt and S^/Nb are almost independent 
of b. This implies that the spin vectors m are hardly influenced by whether the surface is 
smooth or collapsed even though rrij is constrained to be parallel to the triangle i. 

Now, let us turn to the collapsing transition. We show how the IDOF influences the 
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FIG. 4: (a) The internal energy S^/Nb versus 6, and (b) the magnetization M/N^ versus b, which 
were obtained at a = 2, a = l, and a = 0.5. The model is the one defined by Eq.([T]). The range of 
assumed b is spanning the collapsing transition point in each a. 
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FIG. 5: The internal energy S^/Nb versus b at (a) a = 2, (b) a = l, and (c) q = 0.5. 
collapsing transition. The mean square size is defined by 



(7) 



where X is the center of the surface, and is plotted in Figs. E^a)-(c) against b. X"^ in the 
figure was obtained at a = 2, q; = 1, and a = 0.5. We find that the transition point 6c, where 
rapidly varies, moves right on the h axis with decreasing a. In the limit of a = 0, the 
transition point should be 6c 6°, where 6|?^0.77 9] is the transition point of the model 
without the IDOF. This implies that the surface is softened by the interaction between the 
surface and the IDOF, because the stiffness of the surface is considered to be reduced if the 
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transition point be decreases even though the parameter b is itself not always identical to 
the macroscopic bending rigidity of the surface. 
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FIG. 6: The specific heat Cs^ versus b at (a) a = 2, (b) a = l, and (c) a = 0.5. Anomalous peaks 
indicate a collapsing transition between the smooth spherical phase and a collapsed phase. 



The fluctuation of X is defined by 



(8) 



which is expected to reflect how large the surface size fluctuates. Figures [6](a)-(c) show Cx^ 
against 6 at a = 2, a = l, and a = 0.5. As expected from Figs|5](a)-(c), the peak position of 
moves to the right on the b axis as a decreases from a = 2 to a = 0.5. We find also that 
the peak value C™2^ itself increases with decreasing a. This second observation implies that 
the shape fluctuation is slightly suppressed by the interaction between the surface and the 
IDOF, because C^2^ decreases with increasing a. 

In order to see the order of the collapsing transition, we show log-log plots of C^2^ against 
in Figs. [7](a)-(c). The straight lines are drawn by fitting the data to 



y^max 



(9) 



where u is a critical exponent of the collapsing transition. Thus we have 

z/„=2 = 0.75 ± 0.08 (a = 2), 
z/„=i = 1.09 ± 0.05 (a = l), 
1/^^0 5 = 1.13 ±0.06 (a = 0.5). 



(10) 



From the finite-size scaling (FSS) theory, the second and the third results in Eq. fllOl) indicate 
that the transition is of first-order, because u are compatible to z/= 1. On the contrary, the 
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FIG. 7: Log-log plots of the peak value C^^^ against N at (a) a = 2, (b) a = 1, and (c) a = 0.5. 
The straight lines were drawn by fitting the data to Eq.dl]). 
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FIG. 8: The variation of against MCS at (a)-(c) a = 2, (d)-(f) a = l, and (g)-(i) a = 0.5. The 
data were obtained at the transition point be, which depends on both N and a. 

The surface shape is expected to be discontinuously changed at the transition point. 
Then, it is natural to ask what is the Hausdorff dimension H of the surface at the smooth 
phase and at the collapsed phase. We show in Figsl8](a)-(i) the variation of X against MCS 
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obtained at the transition point be with a = 2, a = 1, and a = 0.5. It is almost clear that 
disnontiniioiislv nhanees between the smooth nhase and the nollansed nhase in almost 
all Cc □ 




FIG. 9: Normalized histograms h{X'^) for the distribution of X"^ obtained at (a)-(c) a = 2, (d)-(f) 
a = l, and (g)-(i) a = 0.5. Those h{X'^) correspond to the variations of shown in Figsi8]^a)"(i). 
The dashed vertical lines denote X^f°\ X'i^ and Xll^° , ^mlx°, which will be shown in Tabled 

Figures [9](a)-(i) show the normalized histogram h{X'^) of variation X^, which was shown 
in FigslH](a)-(i). A double peak structure is clearly seen in h{X'^) at all values of a except 
on the surface of size A^ = 4842 in Figj9t^a). The reason of this is only the size effect; h{X'^) 
has the double peak on sufficiently large sized surfaces. Thus, the collapsing transition is 
confirmed to be first-order even at a = 2, where the first-order transition was not confirmed 
by the FSS analysis for C^f^ in Eq.^. 

The double peak in h{X'^) allows us to calculate the mean value X'^^^° of in the 
smooth phase and the mean value '^"^ in the collapsed phase. These mean values can be 
calculated by assuming the lower (upper) bound X^j^ (X^^^^ ) and by averaging X^ in the 
range < X^ < Xl^'^^™°\ which includes one of the two peaks. The dashed lines 

drawn vertically in Figs|9](a)-(i) denote the lower and upper bounds; four dashed lines in each 
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TABLE I: The lower bound X^f°' and the upper bound ^max collapsed state, and the 

lower bound ^mf™° ^^^^ upper bound ^ma™° obtaining the mean value X^(smo) in the 
smooth state. 



a N 
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figure respectively correspond to X^f°\ X^™', and Table U show the lower 

and the upper bounds X^^^, X"^^^ including those shown as dashed lines in Figs|ni^a)-(i). 

Figures [TO](a)-(c) show versus X in a log-log scale, where X^ were obtained by aver- 
aging the data between the lower and the upper bounds listed in Table [H The error bars in 
Figs ITor a)-(c) are the standard deviations. The straight lines were drawn by the power fit 
of the form 

X^ ~ iV2/^, (11) 
where H is the Hausdorff dimension. Then we have 

//smo = 1.85 ± 0.06, Hcoi = 1.90 ± 0.17 (a = 2), 

i/smo = 2.11 ±0.17, //col = 2.37 ±0.43 (a = 1), (12) 
i/smo = 2.20 ± 0.20, H^oi = 2.60 ± 0.52 (a = 0.5). 
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FIG. 10: Log-log plots of the mean values X"^ against obtained in the smooth phase and the 
collapsed phase at (a) q = 2, (b) q = 1, and (c) a = 0.5. The mean values were obtained by 
averaging between the lower bound X^^j^ and the upper bound ^maxi which are shown in Table 
IHand also indicated by vertical dashed lines in Figsl9]||a)-(i). 

The results H^mo in the smooth phase are almost identical to the expected value H = 2, 
which is the topological dimension of surfaces. We see that the results i^coi in the collapsed 
phase obtained at a = 2 and a = 1 are also almost identical to H = 2, and that Hcoi at 
a = 0.5 remains in the physical bound H = 3. In the limit of a ^ 0, iJcoi is expected to be 
if = 2.59(57), which was obtained in the case of a = in [9]. The result ifcoi at a = 0.5 in 
Eg. (1121) is almost identical to this value. Not only the smooth phase but also the collapsed 
phase is therefore considered to be unaffected in the presence of the IDOF in the range of 
small to medium a. 

C. Surface fluctuations 

The bending energy 5*2 can reflect how smooth the surface is, and therefore we plot S2/NB 
versus b in Figs inT a)-(c). where Nb is the total number of bonds. The results in Fig JTTT a) 
at a = 2 are independent of and hence indicate that the surface fluctuation is suppressed 
and the phase transition disappears. On the contrary, S2/NB shown in Figs ITlT a) and [TlT b) 
varies rapidly with increasing A^, and this is considered to be a signal of phase transition 
although a discontinuity can not be seen in those S2/Nb- 

In order to see the order of the transition more clearly, we plot in FigsJTWa)-(c) the 
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specific heat for tlie bending energy 5*2, wliicli is defined by 

Cs,= ^^{{S2-{S2)f). (13) 

The anomalous behavior seen in indicates a phase transition between the smooth phase 
and the collapsed phase. However, as we see in Fig JTW a). the peak value C^^^ remains 
constant even when increases. This is consistent to the behavior of S2/NB in FigJTlTa). 

The peak value Cg"^"^ of is plotted against in a log- log scale in FigsJT3ra)-(c). It is 
apparent that C^^^ in FigIT3l(a) stops growing with increasing A^, which is in sharp contrast 
to the cases in FigsIT3l(b).(c). The straight lines in FigsJT3](b) and fTST c) are drawn by fitting 
the data to the form 

C^^- ^ AT-, (14) 

where a is a critical exponent of the transition. The largest three data are used in the fitting 
in the case of a = 0.5 in FigsJT3](c). Thus, we have 

(7 = 0.63 ±0.04 (a = l), 
(7 = 0.84 ±0.07 (a = 0.5). (15) 

The first result in Eq. ffTSj) implies that the transition is of second order (or continuous) at 
a = 1, because a is considered to be cr < 1. From the finite-size scaling (FSS) theory, we know 
that a <1 corresponds a continuous transition. On the contrary, the transition appear to 
be discontinuous at a = 0.5, because the second result in Eq. (|T5l) is considered to be almost 



S2/NB S2/NB S2/NB 




0.66 0.68 0.7 ■ 0.74 0.75 ' 0.76 0.77 



(a) b (b) b (c) b 

FIG. 11: The bending energy S2/NB versus h at (a) a = 2, (b) a = l, and (c) q = 0.5, where Nb is 
the total number of bonds. 
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a=0.5 



N= 14442 
N=8412 




0.66 0.68 0.7 0.72 0.74 0.75 
(a) ^ (b) b 



0.76 0.77 
(c) b 



FIG. 12: The specific heat €$2 versus h at (a) q = 2, (b) q = 1, and (c) a = 0.5. Anomalous peaks 
(-vmax indicate the existence of phase transition for surface fluctuations. 




FIG. 13: Log-log plots of Cg^^ against N at (a) a = 2, (b) a = l, and (c) a = 0.5. The straight lines 
in (b) and (c) were drawn by fitting the data (the largest three data in (c)) to the form of Eq.([T^ 



equal to a = l. We expect that the order of transition turns to be discontinuous when a is 
reduced, because a discontinuous transition can be seen in the case of a = 9|]. 

Finally, we show Si/N versus h in FigsJTH(a)-(c). Si/N is expected to be Si/N c^3/2 
because of the scale invariant property of the partition function. Consequently, we can use 
this relation to see whether the simulations are performed successfully. It is easy to see that 
the relation is satisfied in almost all cases except in the smooth phase at a = 2; however, the 
deviation is still very small compared to the value Si/N. Therefore, we consider that the 
simulations are performed correctly. 
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1.49 

I 
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0:N= 14442 
A:N=8412 
□ :N=4842 
V :N=2562 



Si/N 

1.51 



1.5 



1.49 



0.66 0.68 0.7 



0:N=14442 
A:N=8412 
□ :N=4842 
V:N=2562 




. 0.74 . 0.75 0.76 . 0.77 

(b) (c) 



FIG. 14: The Gaussian bond potential Si/N versus b at (a) a = 2, (b) a = l, and (c) a = 0.5. The 
expected relation S'i/A^~3/2 is satisfied. 

V. SUMMARY AND CONCLUSION 



We have investigated an interaction between the tilt order and the shape of surfaces of the 
conventional surface model of Helfrich, Polyakov and Kleinert by Monte Carlo simulations 
on triangulated spherical lattices. The purpose of this study is to see how the tilt order 
influences the collapsing transition and the surface fluctuation transition, both of which 
were reported to be of first-order in the conventional surface model jsl. The Hamiltonian 
of the model in this paper is defined by a linear combination of the HPK Hamiltonian and 
that of the XY model. The unit vector nij of the XY model is defined on the triangle i by 
projecting a three-dimensional vector on the triangle, where the three-dimensional vector 
is assumed to represent a lipid molecule usually called a director. Since the vector nij is 
parallel to the triangle i, the XY model is not identical to the naive XY model defined on 
the planar lattices. The parameter a, which is the coefficient of the XY Hamiltonian, is 
assumed to a = 2, a = 1, and a = 0.5 in the simulations. In the case of a = 2 the vectors m 
are considered to be in a relatively ordered state, while in the cases a = l and a = 0.5 they 
are in disordered states. 

We find that the variables mj change depending only on a and are almost independent of 
whether the surface is smooth or not. In fact, the internal energy S-^/Nb and the magneti- 
zation M/Nt remain almost constant in the range of h including the collapsing transition 
point. 

It is also observed that the collapsing transition is not so strongly influenced by the tilt 
order. The transition is slightly softened in the presence of tilt order, however, it remains 
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in first-order and occurs almost independent of a at least up to a = 2. Furthermore, the 
collapsing phase at the transition point is characterized by a physical Hausdorff dimension, 
i.e., ifcoi<3 in all cases a = 2, a = l, and a = 0.5. This result is consistent with the physical 
Hausdorff dimension at a; = reported in Q]. 

On the other hand, the transition of surface fluctuations is influenced by the tilt order. 
The transition appears to remain discontinuous at a = 0.5, where the variables rrij weakly 
correlate with each other and are relatively at random. As the coefficient increases to 
a = l, the transition is softened and turns to be a continuous one. Moreover, the transition 
disappears and turns to be a higher-order one as a increases to a = 2. This result leads us to 
conclude that the collapsing transition is not always accompanied by the surface fluctuation 
transition in the surface model with internal degrees of freedom such as the tilt order. 

Finally, we comment on whether the phase structure is influenced by the singular points 
of m, which appear even at zero temperature due to the surface topology. The exponents 
V shown in Eq.Q and the Hausdorff dimension H in Eq.( fT2l) . both of which characterize 
the collapsing transition, are not influenced by the singular points. In fact, the collaps- 
ing transition is almost independent of the variables m as we have shown in this paper. 
However, it remains unclarified whether or not the exponent a in Eg. (1151) is influenced by 
the singularity; more precisely, it remains to be studied whether or not the softening of 
the surface fluctuation transition is caused only by the singularity of m or the topology of 
surface, although KT transition itself is expected to be independent of the surface topology 
as discussed in Section HV A[ 

This work is supported in part by a Grant-in-Aid for Scientific Research from Japan 
Society for the Promotion of Science. 
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